Quantum Dot between Two Superconductors 
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Novel effects emerge from an interplay between multiple Andreev reflections and Coulomb in- 
teraction in quantum dot coupled to superconducting leads and subject to a finite potential bias 
V . Combining an intuitive physical picture with rigorous path integral formalism we evaluate the 
current / through the dot and find that the interaction shifts the subharmonic pattern of the I — V 
curve is shifted toward higher V. For sufficiently strong interaction the subgap current (at eV < 2A) 
is virtually suppressed. 
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Recent progress in nanotechnology enables the fabri- 
cation and experimental investigation of superconduct- 
ing contacts of atomic size with few conducting channels 
[Q,^ . Transport properties of such systems are essentially 
determined by the mechanism of multiple Andreev reflec- 
tions |3| (MAR) which is responsible for Josephson cur- 
rent as well as for dissipative currents at subgap voltages. 
Theoretical analysis of MAR and current-voltage charac- 
teristics in small superconducting junctions is reported 
in a number of papers |^-^. In these works, an essen- 
tial ingredient is the assumption that electron-electron 
interaction inside the contact can be neglected. It might 
indeed be justified provided a metallic contact is suffi- 
ciently large and/or strongly coupled to massive super- 
conducting leads. 

However, in very small contacts (quantum dots), the 
Coulomb interaction is not effectively screened, hence it 
is expected to substantially affect transport properties 
of the system. For instance, it is well known both from 
theory and experiment that tunneling through a 
quantum dot between superconductors can virtually be 
suppressed due to Coulomb effects. Thus, to the fasci- 
nating physics of SNS and SIS junctions, one should add 
that of an SAS junction composed of superconducting 
leads coupled by an interacting quantum dot. 

In the present work, the physics of interplay between 
MAR and interaction effects in SAS junctions subject to 
a finite bias is exposed. It encodes the salient features 
of superconductivity, strong correlations and non-linear 
response. A simple intuitive physical picture is combined 
with a rigorous path integral technique by which irrele- 
vant degrees of freedom are eliminated and an effective 
action is constructed (in the spirit of Feynman- Vernon in- 
fiuence functional pf). Similar ideas proved to be useful 
elsewhere, see e.g. ]10| , pT| . In the present context they 
have been applied for the relatively simple case of an 
SAS junction at zero bias, focusing on the equilibrium 
Josephson current ||l^,|l^. Our main achievements are: 
a) Derivation of a tractable expression for the non-linear 
tunneling current in the presence of interactions, b) pre- 
diction of novel physical effects pertaining to the I — V 
curve of an SAS junction at sub-gap bias. 

Let us then commence with a simple and physically 



transparent picture of an interplay between MAR and 
Coulomb effects. Consider a quasiparticle (hole) which 
suffers n Andreev reflections inside the superconducting 
junction, thereby gaining an energy neV, where V is the 
voltage bias. As soon as nev = 2A the quasiparticle 
leaves the junction and does not contribute anymore to 
the subgap current. Hence, the number of Andreev re- 
flections n for a given voltage is n ~ 2A/eV. 

Assume now that the Coulomb interaction inside the 
junction is switched on. For our qualitative discussion it 
suffices to account for it in terms of an effective capac- 
itance C and its related charging energy Ec — e^/2C 
At T = and for eV < Ec a single electron tunneling 
(and, hence, also MAR) is blocked, so in what follows 
we will consider the case eV > Ec > 0. In order to 
leave the junction the quasiparticle should gain an energy 
neV equal to 2A + {n + l)Ec- The last term originates 
from the fact that during the MAR cycle with a given n 
the charge {n+ l)e is transferred between the electrodes. 
Hence, an additional energy (n -I- l)Ec should be paid. 
The above condition immediately fixes the number n at 
a given voltage: 



2A + Ec 



eV -Ec 



(1) 



Thus, in the presence of Coulomb interaction quasiparti- 
cles spend more time inside the junction and suffer more 
Andreev refiections. Another obvious observation is that 
at low temperature T, the transfer of the charge (n-f l)e 
is blocked by interaction at voltages eV < (n + \)Ec- 
Combining this observation with (|l|) one arrives at the 
condition 



zV<eVth^Ec I l + \/l + ^j 



(2) 



under which the MAR current is suppressed due to 
Coulomb repulsion. For Ec <C A the voltage threshold 
is eVth ^ \J2l\Ec > Ec, i.e. in this case MAR should 
be blocked even at voltages much higher than Ec / e- For 
eVth > 2A eq. (|) yields 

Ec > 2A/3. 
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For such values of Ec one expects the subgap current to 
be fully suppressed due to Coulomb interaction. 

Let us now recall that the subharmonic peaks occur 
on the I — V curves each time the MAR cycle with a 
given n becomes impossible. Without interaction these 
peaks are located at voltages Vn = 2 A/ en. It follows 
immediately from the above discussion that in the pres- 
ence of Coulomb interaction the peaks should be shifted 
to higher voltages. From eq. (|]) one finds 



^ ^ 2A + Ec 

e en 



(3) 



i.e. one expects the subharmonic peaks to be shifted by 
SVn = Ec{l + towards larger V as compared to the 
noninteracting case. 

Thus, already a naive analysis of the interplay between 
MAR and Coulomb effects allows one to predict several 
novel effects which can be experimentally tested. To put 
these qualitative arguments on a firm basis we formulate 
below a realistic model of an SAS junction, and proceed 
with a rigorous calculation of the I — V characteristics. 

The model and basic formalism. Consider, in two di- 
mensions, a quantum dot at r = weakly coupled to (half 
planar) superconducting electrodes. The Hamiltonian of 
the system is decomposed as. 



H = Hi 



Hi 



H, 



dot 



Ht 



(4) 



The Hamiltonians of the left (x < 0) and right (x > 0) 
superconducting electrodes have the standard BCS form 



dr[^l{r)a^)^Ar) 



AvI/Ur)vI/t (r)vI/,^(r)*,T(r)]. 



(5) 



Here ^f^^ {"^ja) are the electron creation (annihilation) 
operators, ^(V) = — V^/2m — /i, and j = L,R for left 
and right electrodes. The dot itself is modeled as an 
Anderson impurity center with Hamiltonian 



(6) 



where Cj. and are dot electron operators. The im- 
purity site energy eo (counted from the Fermi energy /i) 
is assumed to be far below the Fermi level eo < 0. The 
presence of a strong Coulomb repulsion U > —eo between 
electrons in the same orbital guarantees that the dot is 
at most singly occupied. 

Electron tunneling through the dot is described by the 
term. 



j=L,R <T 

where T is an effective transfer amplitude. 



(7) 



The dynamics of the system is completely contained 
within the evolution operator on the Keldysh contour 
K jlj] which consists of forward and backward oriented 
time paths. Its kernel J is given by a path integral. 



J= V^V^VCVCe^piiS), 



(8) 



over Grassman fields corresponding to the fermion opera- 
tors, with ^ = (^f^^^, ^l^j^, vp^jp ^f^^) with obvious defini- 
tions for 5*, C and C. Moreover, S — Ldt is the action 
and L is the Lagrangian pertaining to the Hamiltonian 



In order to avoid dealing with fields defined on both 
branches of the Keldysh contour one performs a rotation 
C ^ c and ^ ip in Keldysh space: 



= Co 



c^QC; Q^^ 



1 -1 
V^l 1 1 



(9) 



and similarly for ■0 and ip. Here ct^ is the third Pauli 
matrix operating in Keldysh space. The new Grassman 
variables c, c, ip, ip are now defined solely on the forward 
time branch. Averages of the corresponding products of 
these fields determine the standard 2x2 Green-Keldysh 
matrix composed of retarded (G^), advanced (G^) 
and Keldysh {G^) Green functions which, in turn, are 
2x2 matrices in spin (Nambu) space. 

The path integral (||) is expressed in terms of the new 
Grassman variables in the same way, and the action S is 
now defined as S* = S'dot + Sq[iJj,iI>], where 



S'dot = I dt 



.d_ 

dt 



U 



[ccf 



(10) 



so^ [dtY: If 



drV;, (r,i)G7V, (r, t) 

j=L,R '--J 

+ {Ti:j{0,t)T,c{t)+c.c. 



(11) 



where e = eo -I- U/2 and the Pauli matrices r^^y^z act in 
Nambu space. The operator G^^^ has the standard form 

where t± = [t^ ±iTy)/2 and A^^n are the (spatially con- 
stant) BCS order parameters of the electrodes. 

Effective action and transport current. The basic al- 
gorithm of our approach is to integrate out the electron 
variables in the superconducting electrodes which play 
the role of an effective environment for the dot. This 
procedure yields the influence functional F[c.^c] for the 
c- fields in the dot: 

F = exp(iS'e„v[c,c]) = / I?V;i?Vexp(zS'o[Vi,V']), (13) 
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which is evaluated exactly. Gaussian integration in jl^ ) 
is carried out separately for L- and i?-electrodes. 

Let us consider, say, the left superconductor and omit 
the subscript j = L for the moment. The first step is 
to integrate out the fermion fields inside the supercon- 
ductor thereby arriving at an intermediate effective ac- 
tion in terms of the fermion fields defined on the surface 
a; = 0. It is useful at this point to Fourier transform the 
fields ip^x, y) along the (translationally invariant) y di- 
rection. The problem then reduces to a one dimensional 
one with fermion fields ipkix) where k is the quasiparti- 
cle momentum in the direction normal to x. In order to 
evaluate the Gaussian integral we will look for a saddle 
point field ipkix) defined by G~^{£,x)ipk{x) = 0, where 

= -{l/2m){d'^/dx'^) — and fi^ = fi ~ k"^ /2m. 

Decomposing -0 into bulk and surface fields ipkix) = 
tp\ {x) -\- ipk (0) and integrating out V'fc {x) we arrive at the 
intermediate effective action 5 of a superconductor lead 
expressed only via the T/j-fields at the surface, 

S^iJdtJ dt'J2Y'^k{0,t)T,g{t,t')M0,t'). (14) 



Here Vx = \/2/ife/m and 

git,t') = j 5(6)e""(*-* )^e-^^, (15) 

is the Green-Keldysh-Eilenberger matrix of the (left) su- 
perconducting electrode 
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(16) 



(p{t) = ipo + 2e J V{ti)dti is the time-dependent phase 
of the superconducting order parameter and V(t) is the 
electric potential of the electrode. The Fourier trans- 
formed retarded and advanced Eilenberger functions have 
the standard form 



(e±zO)T^ + z|A|r^ 
V(e±*0)2-|A|2' 



(17) 



and g^ = (g-" - g"^) tanh(e/2T) is the Keldysh function. 

The second step in our derivation amounts to integrat- 
ing out the ^-fields on the surface. The integral 



r»V;(0)X>V(0) exp(iS' + i dt{TiPkiO)TzC + c.c.)) (18) 



can easily be evaluated. Carrying out exactly the same 
procedure for the right electrode, making use of the iden- 
tity g\ p = ^ and adding up the results we obtain 

5e„v ^iT I dt ( dt'c{t)T,g+{t,t')c{t'). (19) 



Here and below we define F = 4^j,T^/u2: and g± 
igL±gB.)/2. 



Eq. ( p^ ) is one of our central results. It enables the 
expression of the kernel J (||) solely in terms of the fields 
c and c: 



J = / 'Dc'Dcexp{iS dot +iS'onv), 



(20) 



where S'dot + 5'oiw = '5'cff [c, c] represents the effective ac- 
tion for a quantum dot between two superconductors. 

In order to complete our derivation let us express the 
current through the dot in terms of the correlation func- 
tion for the variables c and c. Starting from the general 
expression for the current and representing the correlator 
for the ip-Helds in terms of that for the c-fields we find. 



Y Tr [g( cc) I A- 



h.c.]. 



(21) 



Thus, the problem of calculating the current through an 
interacting quantum dot is reduced to that of finding 
the correlator (cc) in the model defined by the effective 
action S^tf = S'dot + •S'onv (|lO|), (|T^). It should be em- 
phasized that our approach is appropriate for studying 
both equilibrium and nonequilibrium electron transport. 
In the noninteracting limit [/ — s- the results of previous 
studies can be easily recovered within our formalism. 

Mean field approximation. Consider now the case 
U ^ and decouple the interacting term in ( [lO|) by 
means of Hubbard-Stratonovich transformation [|l2|jl3[ 
introducing additional scalar fields 7± . The kernel J now 
reads. 



J = J VcVcV-f+Vj^ exp(zS'[7] + iS^slu^o) 



dt 



c"f+axC + cj-c - 



(22) 



(23) 



Here we will assume that the effective Kondo temper- 
ature [0 Tk —^/UT exp [— 7r|eo|/2F] is smaller than the 
superconducting gap A. In this case interactions can 
be accounted for within the mean field approximation. 
The fields 7± in ( p3|) are considered as time-independent 
parameters determined self-consistently from the saddle 
point conditions 5J/Sj± — 0: 



7+ 



U 



dt < cc>, 7_ 



U 
2" 



dt < coxc > . (24) 



As it turned out from our numerical analysis the effect 
of the parameter 7-|_ is merely the renormalization of the 
tunneling rate F. Absorbing 7^ -terms in F we arrive at 
the final effective action of our model 



^cff[7]= / ^ I de'cM{e,e')c, 



(25) 



M{e, e') = S{e - e'){e + 7- - T,~e) + iT,Tg+{e, e'). (26) 
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FIG. 1. The I-V curves for an SAS junction at subgap volt- 
ages. The parameters are eo = —1.5, T = 0.6, U = 2.4 (sohd 
curve) and 2.7 (dashed curve). 

Subgap current in SAS junctions. In order to find the 
correlator (cc) = iM~^ and the current (|^) we numeri- 
cally inverted the matrix ( p6| ) and simultaneously solved 
the self-consistency equation for 7_ (|2^). The resulting 
I — V characteristics for an SAS junction in the pres- 
ence of Coulomb interaction are displayed in Fig. 1 . One 
observes all the main features predicted within our sim- 
ple picture of an interplay between MAR and Coulomb 
interaction: (i) at relatively low voltages V < Vth MAR 
current is essentially suppressed due to interaction, (ii) 
for higher voltages (but still eV < 2A) MAR is possible 
and results in a nonzero subgap current which increases 
with V and (iii) the subharmonic peaks in the differen- 
tial conductance occur and are shifted to higher voltages 
as compared to the noninteracting case. An increase of 
U results in a stronger current suppression and a more 
pronounced shift of the subharmonic peaks. Close to the 
gap edge eV — 2A the current shoots up sharply. 

The parameters used in our numerical analysis are cho- 
sen in a way to observe all the key features (i), (ii) and 
(iii). It is interesting to quantitatively compare the re- 
sults presented in Fig. 1 with the predictions (0)-(H) of 
an oversimplified "iJc-based" model. Let us estimate 
the effective value of Ec (which is, strictly speaking, a 
function of U, ep T and V in our calculation) with the 
aid of eq. (||) and the I ~ V curves of Fig. 1. We 
find Ec » 0.2A for U = 2.4A and Ec « 0.25A for 
U = 2.7A. Obviously, these values of Ec are smaller 
than 2 A/3 and, hence, a finite subgap conductance is 
expected at eV ^ \/2AEc ^ Ec ■ This is precisely what 
we observe in Fig. 1. For such values of V and Ec eq. 
(|l|) yields n < 3, i.e. only two subharmonic peaks (with 
n = 2, 3) can occur. This is exactly the case in Fig. 1. Fi- 



nally, substituting the above values of Ec into eq. (Q) we 
can estimate the magnitudes of the peak shifts 6Vn- For 
U = 2.4A we find SVn ~ 0.3A for n = 2 and - 0.26A for 
n = 3. Analogous values for U = 2.7 A are respectively 
~ 0.38 A and ~ 0.33 A. These values are in a reasonably 
good agreement with our numerical results. 

In summary, we presented a detailed analysis of an 
SAS junction at finite bias and derived its effective action 
using Keldysh path-integral techniques. Our approach 
applies for both equilibrium and nonequilibrium current 
transport in the presence of interactions. The repulsive 
Coulomb interaction leads to novel effects in the pattern 
of the subgap current. In particular, it shifts the peaks 
of the differential conductance toward larger bias. When 
the interaction is sufficiently strong the subgap current is 
highly suppressed. Our theoretical predictions can be di- 
rectly tested in experiments with superconducting quan- 
tum dots. 
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